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One of the major themes in the physics of condensed matter is unconventional superconduc¬ 
tivity in correlated electronic systems. The symmetry structure of Cooper pairs is thought to 
be the key to the understanding of the pairing mechanism of their superconductivity 12 . In 
iron pnictide superconductors, the precise symmetry of the order parameter as well its evolu¬ 
tion with doping remains highly controversial 3 ! There is no general consensus on the nature 
of pairing in iron-based superconductors leaving the perspectives ranging from s ++ wave, to 
s ± , and d wave. Here we show that doping induced a sudden change of the pairing symmetry 
from the s ++ - to d -wave in the BaFe 2 -a;Ni x As 2 system with an ^ state near the crossover 
boundary at x — 0.15. The presence of the s ± state is of particular scientific interest because 
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an intermediate phase breaks the pairing symmetrj^" 9 when the s- and d- wave phases come 
together to interfere. For x < 0.15, the temperature dependence of the London penetration 
depth, lower critical field, and specific heat in the superconducting (SC) state provides strong 
evidence for a two-band s-wave order parameter. Upon doping for x > 0.15, we find the tem¬ 
perature and magnetic-field contributions to the specific heat in T 2 and \/~R, respectively, 
indicating line nodes in the SC gap. These findings thus clarify the symmetry nature as well 
as the evolution of the SC order parameter with doping in pnictide superconductors. 

Conventional phonon-mediated superconductors and unconventional cuprate superconductor^ 
are well characterized by distinct s-wave and d -wave pairing symmetries with nodeless and nodal 
gap distributions, respectively. In conventional superconductors, the electron-phonon interaction 
gives rise to the attraction between electrons near the Fermi-surface (FS) with opposite momenta 
and opposite spins, which eventually causes superconductivity. On the other hand, it has gener¬ 
ally been believed that the electron-electron coupling leads to the formation of Cooper pairs in 
cuprates and pnictides 3 . However, iron pnictides present a rich phase diagram wherein supercon¬ 
ductivity coexists and competes with the spin density wave, resulting in unconventional pairing 
mechanism^ 10 !. Although there is a general consensus that spin fluctuations play an important role 
in the formation of Cooper pairs, many aspects such as the role of magnetism, the nature of chem¬ 
ical tuning, and the resultant pairing symmetry remain unsettled 2 3 . In a d -wave pairing symmetry, 
the order parameter changes sign with angles in the x-y plane, forcing the gap to go to zero along 
diagonal directions (k y = ±/c T ) 11 . In the case of s-wave symmetry, the gap in general will not 
have nodes. In the momentum space, an s-wave SC gap has an isotropic magnitude in all direc- 
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tions, and it has a fixed phase in all directions (see the inset of Fig. 1). More specifically, it is 
additionally experimentally shown that the electronic specific heat, angle-resolved photoemission 
spectroscopy (ARPES), and scanning tunneling micrposcopy (STM) studies of Ba(Fci __ x Co x ) 2 As 2 
agree on an isotropic s-wave pairing symmetry 12 14 } while Raman scattering and heat transport 
measurement 15 ^ arc interpreted in terms of nodes. Theoretically, an s =t -wave pairing symmetry 
(nodeless gap function with opposite signs of the order parameters for electron and hole pock¬ 
ets) has been suggested 17 . Upon further doping, electron hole bands disappear, making the sit¬ 
uation in the overdoped regime not so clear, though several forms of the order parameters have 
been proposed such as extended s, d, and s + id -wave pairing symmetries for heavily overdoped 
compounds 18 19 . Therefore, the SC gap structure of the iron-based superconductors is not univer¬ 
sal. It is quite different from that of cuprates in which almost all have a nodal pairing state 11 . Such 
scattered pairing symmetries and various interpretations occur partly due to a sensitive dependence 
on measurement probes and material stoichiometry. 

For clarifying this situation, we chose an electron-doped BaFe 2 - x Ni x As 2 system over the 
whole doping regime from its parent to an overdoped non-SC metal. The specific heat was system¬ 
atically measured for 13 samples as marked by the arrows in Fig. 1. The specific heat experiment 
has probably the best energy resolution of all experimental probes for distinguishing nodes in the 
gap. Near the absolute zero point, the bulk thermodynamic quantity is a sensitive probe of low- 
energy excitations of a complex quantum system. Low-energy excitations contain useful informa¬ 
tion about the nature of the ground state and can help to evaluate the presence or absence of nodes 
in the SC order parameter. For an s-wave superconductor, the electronic contribution to the specific 
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heat is given by: C e i oc exp where A 0 is the energy gap of the SC density of states (DOS). For a 
nodal gap symmetry, Sigrist and Ueda 20 showed that the line nodes in the gap are derived from the 
temperature dependence of the specific heat ocT 2 . Here, we present the experimental evidence for 
the crossover of the pairing symmetry from nodeless s-wave to node d -wave of BaFe 2 _ a .Ni x As 2 
around x = 0.15. The results are in good agreement with the early experiments on the optimally 
doped compound 21 22 and the penetration depth measurements on overdoped materials 23 . 

To have full information on physical properties, we traced in Fig. 1 the doping dependence 
of structural, magnetic, and SC properties of BaFe 2 - x Ni x As 2 single crystals as a function of Ni 
content. The undoped parent (x = 0) undergoes both an antiferromagnetic (AFM) ordering and 
a structural transition from paramagnetic tetragonal to an AFM-ordered monoclinic phase at T 0 
= 136(3) K with a first-order characteristic. Upon the Ni substitution, the sharp first-order struc¬ 
tural/magnetic anomaly of the parent compound gradually broadens and shifts to lower temperatures 24 . 
The commensurate static AFM order changes into transversely incommensurate short-range AFM 
order near optimal doping that coexists and competes with superconductivity. Upon further in¬ 
creasing the Ni concentration, the anomaly is no longer visible in the specific heat data, while 
superconductivity maintains at slightly lower temperatures, i.e., T c = 0.85 K and 2.8 K for x = 
0.22, and 0.20, respectively. 

Figure 2a illustrates the electronic specific-heat contribution of superconducting sample with 
x = 0.096, obtained by subtracting the lattice specific heat from the raw data following the proce¬ 
dure given in ref. [24]. The entropy conservation required for a second-order phase transition is 
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fulfilled. It is obvious that the SC transition at T c is well pronounced showing a sharp jump in C e \. 
We observed that superconductivity emerges in the magnetic phase for x <0.10 accompanied by 
increasing T c . This increase is accompanied by a rise of the electronic coefficient 7 ,, of the specific 
heat in the normal-state from 5.1 mJ/mol K 2 at the undoped sample to about 24mJ/mol K 2 to the 
optimal concentration, which we attribute to the closing of the SDW gap with Ni doping 24 . Our 
data stresses that the decrease of with increasing x is linear for x > 0 . 10 , while the residual 
electronic specific heat y r linearly rises. This explains the increase of T c and AC c \/x n T c with Ni 
doping in the underdoped region, implying that a larger fraction of the FS becomes available for 
superconductivity. However, the almost linear temperature dependence of C e \/T of the SC samples 
indicates that the specific heat data cannot be described by a single BCS gap. In order to illustrate 
this we show a theoretical BCS curve with A = 1.764 k B T c = 2.23 meV in Fig. 2a. A precise 
description of the experimental data for the two-gaps s-wave model is obtained by using values of 
Ai(0) = 1.8k s T c and A 2 (0) =0.74k s T c . 

The magnetic field dependence of the specific heat through a vortex excitation in a mixed 
state is another excellent, sensitive independent testing method for identifying the gap structure of 
superconductors. It has been well demonstrated that 7 (H) oc H for the isotropic s-wave because 
the specific heat in the vortex state is dominated by the contribution from the localized quasiparti¬ 
cle in the vortex core 25 . It is clear from Fig. 3 that the magnetic field enhances the low-temperature 
specific heat continuously, indicating the increase of the quasiparticle DOS at the FS induced by 
a magnetic field. The fact that the low-temperature specific heat data exhibits a linear behavior 
with magnetic field at low temperatures without any upturn indicates the absence of Schottky-like 
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contributions in our samples. The data in Fig. 3a roughly increases linearly with the applied mag¬ 
netic field. Therefore, the observations of the temperature and field dependencies of the specific 
heat clearly show that the SC gap in the underdoped compounds is nodeless. The temperature 
dependence of the London penetration depth in the underdoped compounds provides strong evi¬ 
dence for a two-band s-wave order parameter 24 . In addition, the lower critical field H c \ values in 
BaFe 2 _ x Ni x As 2 samples are measured by following the procedure described in ref. [24]. The cor¬ 
rected values of H c i are presented and associated with two-gap features. The s ++ wave is therefore 
identified to be the symmetry structure of the underdoped compounds. Additionally, the two-gap 
model also reproduces the electronic contribution to the specific heat of other doped samples well. 
The similar s ++ pairing symmetry has been observed in a phonon-mediated superconductor MgB 2 
which has two s-wave energy gaps with the same phase in the cr and 7 r bands 26 . 

However, the energy spectrum of quasiparticles E k in a singlet SC state can be given by 
E k = + A|, where £& represents the energy of an electron relative to the chemical potential 

with momentum k. On the other hand, unconventional superconductivity is mostly characterized 
by the anisotropic SC gap function with nodes along certain directions in the momentum space. In 
the case of gap nodes, quasiparticles are generated to the largest extent in the vicinity of gap nodes. 
Moreover, the quasiparticle excitation spectrum takes the form E k = Jl ^/upA; 2 + r'^/c 2 , where k± 
and k\\ are wavevectors perpendicular and parallel to the FSs, respectively, vf is the renormalized 
in-plane Fermi velocity at the position of the node, while va ~ <9 A / hdk is the slope of the gap at 
the node associated with the dispersion of the quasiparticles along the FS 27 . This leads to a linear 
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dependence of the DOS with energy. For the two-dimensional case it reads: 


where i denotes the sum over all nodes. The observed T 2 -behavior in Fig. 2c suggests that this 
condition is fulfilled for the overdoped samples i.e., x = 0.18, 0.20, and 0.22. Using Eq. 0 we 
also assumed that the energy gap at zero temperature, A 0 , is equal for all bands, while is a 
renormalized normal-state Sommerfeld coefficient, being an average value over all FSs. 

The SC order parameter C/T oc H 05 holds, if the energy associated with the Doppler shift 
of d- wave quasiparticles in the vicinity of vortices E H = hoyi"II/ To) 0,5 is comparable to E r = 
kbT , where Tq is the flux quantum 27 . According to our data shown in Fig. 3c for x = 0.18, the 
field dependence of C/T follows a H 0 - 5 - law which corresponds to the theoretical prediction 27 . 
Therefore, the observed T 2 term together with the observed H 0 5 behavior of the specific heat in 
the SC state for the overdoped samples evidences d- wave superconductivity in almost all FS sheets. 
In that case the gap would necessarily have four line nodes that run vertically along the c axis. In 
such a nodal structure, zero-energy nodal quasiparticles will conduct heat not only in the plane but 
also along the c axis by an amount proportional to the c-axis dispersion of the FS 11 . However, the 
nonlinear 7 (H) found in Fig. 3c cannot be simply attributed to the multi-band effect as observed 
in Fig. 3a in which 7 ( II) is linear. Even if the multi-bands exist, the zero-field specific heat data 
shows a clear flattening in the low-temperature range corresponding to the weak excitation of 
quasiparticles for an s-wave superconductor. This is completely absent for the three compounds of 
x = 0.18, 0.20, and 0.22 (Fig. 2c). 
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In our investigated system, near the boundary between s ++ and d states, we found an 
state. Theoretical 1\ 2K , the ^-pairing state with different gap sizes (A sma ii ^ A i arge ) represents a 
strong field dependence 'y(H) oc \/~H — H. In addition, the most important virtual excitations that 
trigger the ^-pairing are the AFM spin fluctuations. Although .s ± -wave and d- wave are due to the 
Doppler shift effect by the supercurrent circulating around the vortices, the linear field dependence 
in the s ± -wave model is in contrast to the \[U dependence 28 . In Figs. 2b and 3b, we overview 
theoretical calculations for both iA and d- wave approaches for x = 0.15. It is obvious from the 
field dependence data that the s ± -wave model fits well with the field dependence which is taken 
as a signature for a nodal gap structure. In Fig. 2b, we show the temperature dependence of the 
electronic contribution to the specific heat for x = 0.15, obtained by subtracting the lattice specific 
heat as presented in ref. [24]. From an entropy-conserving construction a T c = 13.12 K was found. 
The specific heat jump at T c is 14.5 mJ/mol K 2 at zero field. The jump relative to the normal 
state specific heat is 0.83, which is much smaller than the BCS value 1.43. This indicates that 
the quasiparticles participating in the SC condensation experience strong elastic scattering because 
inelastic scattering usually enhances the jump ratio 28 . However, the value of the x r /x n in the SC 
state for x = 0.15 accounts for 26.6%. This value also supports the presence of strong scatterers 
because weak scatterers are not capable of inducing states inside the SC gap of unconventional 
superconductors 28 . However, in the case of s 1 * 1 superconductivity, 7 r can be understood as arising 
from interband scattering. This could be induced by disorder due to Ni doping, which is pair¬ 
breaking for a sign-reversing order parameter 17 . 

A theoretical BCS curve with A = 1.764kBT c = 2.23 meV is shown in Fig. 2b with a 
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small fraction 7 res /7n = 0.09 of normal electrons. An .‘A-wave pairing model was applied by 
combining the quasiclassical formalism with the first-principles calculation 29 . Moreover, we used 
the d -wave approach to represent the generic behavior of a nodal SC gap state. More obviously, the 
temperature dependence presents a slightly concave-down behavior. This feature is characteristic 
of a typical multigap superconductor, in which a rather small SC gap coexists with large main 
gaps. This behavior suggests not only the existence of a small gap but also that of a medium-gap 
band 29 . In Fig. 2b, the .A approach is close to the weak-coupling single-band BCS result. One of 
the reasons is that the weighting of the small-gap band is small compared to the total one, i.e., the 
ratio N 3 lN t is 0.2, where N t is the total DOS 29 . In addition, the gap-amplitude difference between 
the minimum |A 3 | and the maximum |A 2 | is not so large, where |A 3 | /|A 2 | ~ 0.5. Moreover, 
|A 2 |/T c is found to be « 2, which is bigger than the single-band value 1.76. Thus, the specific 
heat data from the temperature and magnetic field measurements is in favour of the sign changing 
multiband superconductivity of the s ± -wave state for x - 0.15. 

The reasoning for superconductivity may look quite straightforward and the s 1 * 1 -pairing is 
favored when the antiferromagnetism is suppressed. To get .A superconductivity, one has to invoke 
some mechanism to enhance interpocket interactions 8 . However, the majority of researchers be¬ 
lieve that s ± is the right symmetry, even though the structure of .A SC order parameter has turned 
out to be more complex than originally thought. In strongly electron-doped Fe-based superconduc¬ 
tors, d -wave superconductivity comes following .sA but it emerges as the leading SC instability, for 
which the electron-hole interaction is reduced. The observation of a change of pairing symmetry 
in the same system upon doping would be unprecedented and is another reason why researchers 
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are so excited about Fe-based superconductor^. Several groups have argued that a change from 
s-wave to d -wave would probably proceed through an intermediate state such as s + is or s + id 
symmetry, a complex order parameter that breaks time-reversal symmetry and C 4 lattice rotational 
symmetry 5 9 . In the current work the specific heat data was obtained on the samples with sys¬ 
tematic doping variation by using a single bulk detection technique. Our investigated underdoped 
compounds are definitely consistent with nodeless multiband superconductivity. The large T 2 con¬ 
tribution and the \fJl behavior at low temperatures for the overdoped samples provide evidence for 
the presence of line nodes in the SC gap. For the slightly overdoped level (x = 0.15), we find that 
6^ works well from both T and H dependence, indicating that its SC state is most compatible with 
the multiband s 41 wave. However, at low temperatures the .s :± pairing symmetry could coexist with 
other SC phases, i.e., t -type superconducting orders 5 30 , and the relative phase between the two is 
manifesting explicitly the breaking of the time-reversal symmetry. We thus obtained an interesting 
crossover of the pairing symmetry from s ++ to d -wave by passing an ,s ± for a typical iron-pnictide 
family, which is believed to host unconventional pairing symmetry in these superconductors. 

METHODS SUMMARY 

BaFes-.NUAsa (x = 0, 0.03, 0.065, 0.085, 0.092, 0.096, 0.10, 0.12, 0.15, 0.18, 0.20, 0.22, and 
0.25) single crystals were grown by the FeAs self-flux method, details for the growth process and 
sample characterization were published elsewhere 10 24 . The actual Ni level was determined to be 
80% of the nominal level x through the inductively coupled plasma analysis of the as-grown single 
crystals. Magnetization measurements were performed by using a Quantum Design SC quantum 
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interference magnetometer. The low-T specific heat down to 0.4 K was measured in its Physi¬ 
cal Property Measurement System with the adiabatic thermal relaxation technique. Specific heat 
measurements were performed down to 70 mK by using a heat-pulse technique within a dilution 
refrigerator along H || c up to H = 9 T. 
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Figure 1 | Gap symmetry crossover in electronic phase diagram of BaFe 2 _ a; Ni i As 2 . 

The phase diagram obtained from the magnetic and specific heat data, showing the sup¬ 
pression of the magnetic ( T N ) and structural (T s ) phase transitions with an increasing Ni 
concentration and the appearance of the SC transitions. The arrows indicate the various 
doping levels studied in this work. Superconductivity exists in a dome-like region below 
the transition temperature T c . With increasing x, T c rises to reach a maximal value of 20 K 
at optimal doping. We show that the symmetry of the SC state changes from s-wave to 
d-wave (line nodes) and in between an s ± . The inset illustrates the sketch of the SC gap 
of s-wave vs. d-wave gapP. The top row shows the momentum space representation in 
the regions of the Brillouin Zone having a (+) phase (pink) or (-) phase (blue). The bottom 
row illustrates the SC gap on a large Fermi surface. The s-wave has isotropic gap and 
it does not change sign. The d -wave gap changes sign four times as it goes around the 
Fermi surface and as a result the gap necessarily goes to zero at four nodes. 

Figure 2 | Temperature dependence of the specific heat of BaFe^N^As^ The 

measurements were performed in zero magnetic field and temperature down to T = 
70 mK. The dashed lines in (a) for x=0.096 and (b) for x=0.15 represent the theoretical 
curves based on single-band weak-coupling BCS theory, while the dotted lines illustrate 
the d-wave approximation!^. The solid red line in (a) for x=0.096 indicates the two s-wave 
gap model. The dot dashed line (b) for x = 0.15 represents the s ± -wave modeP. The 
dotted lines in (c) for x = 0.18, 0.20, and 0.22 (two samples) are the linear fits to the 
experimental data at low temperatures to extrapolate the residual linear terms at T = 0. 
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Figure 3 | Magnetic field dependence of the specific heat of BaFe 2 _ a NLAs 2 . The 

magnetic field was applied along the c axis for the specific heat measurements, i.e., H || c. 
The solid line in (a) shows the linear fit to the experimental data for £=0.096, while the 
dotted line in (c) represents a fit to 7 (H) = A(/i 0 H) 0 - 5 for x = 0.18. For £=0.1 5 in (b), the dot 
dashed line is the theoretical fit with the s ± -wave model and the dotted line is the fit with 
the d-wave model, respectively. For an anisotropic gap, the specific heat is deviated from 
the field linear dependence. 
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